Thermal rectification properties of multiple-quantum-dot junctions 
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It is illustrated that semiconductor quantum dots (QDs) embedded into an insulating matrix con- 
nected with metallic electrodes and some vacuum space can lead to significant thermal rectification 
effect. A multilevel Anderson model is used to investigate the thermal rectification properties of the 
multiple-QD junction. The charge and heat currents in the tunneling process are calculated via the 
Keldysh Green's function technique. We show that pronounced thermal rectification and negative 
differential thermal conductance (NDTC) behaviors can be observed for the multiple-QD junction 
with asymmetrical tunneling rates and strong interdot Coulomb interactions. 



Records of thermal rectification date back to 1935 
when Starr discovered that copper oxide/copper junc- 
tions can display a thermal diode behavior. 1 Recently, 
thermal rectification effects have been predicted to oc- 
cur in one dimensional phonon junction systems. 2-5 Such 
a thermal rectification effect is crucial for heat storage. 
Scheibner and coworkers have experimentally observed 
the asymmetrical thermal power of the two-dimensional 
electron gas in QD under high magnetic fields. 6 So far the 
rectification mechanism of a single QD is still ambiguous 
owing to the unclear relation between the thermal power 
and the thermal rectification effect. This inspires us to 
investigate whether the QD junction system can act as 
a thermal rectifier. A useful thermal diode to store so- 
lar heating energy not only requires a high rectification 
efficiency but also high heat flow. The later requires a 
high QD density in the QD junction thermal diode. The 
main goal of this study is to illustrate that the multiple 
QDs embedded into an insulator connected with metallic 
electrodes and with a vacuum layer insert can give rise to 
significant thermal rectification and negative differential 
thermal conductance (NDTC) effects in the nonlinear re- 
sponse regime. We also clarify the relation between the 
thermal power and the rectification effect. 

The proposed isulator/quantum dots/vacuum (IQV) 
double barrier tunnel junction system (as illustrated in 
Fig. 1) can be adequately described by a multi-level An- 
derson model. 7 Here, the vacuum layer serves as a block- 
ing layer for phonon contributions to thermal conduction, 
while allowing electrons to tunnel through. We assume 
that the energy level separation between the ground state 
and the first excited state within each QD is much larger 
than k R T, where T is the temperature of concern. There- 
fore, there are only one energy level for each QD. We have 
ignored the interdot hopping terms due to the high poten- 
tial barrier separating QDs. The key effects included are 
the intradot and interdot Coulomb interactions and the 
coupling between the QDs with the metallic leads. Using 
the Keldysh-Green's function technique, 8 the charge and 
heat currents through the junction can be expressed as 
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l/(exp^~^ L ^^ kBTL ^ + 1) is the Fermi distribution 
function for the left (right) electrode. The chemical po- 
tential difference between these two electrodes is related 
to the bias difference hl — = eAV created by the 
temperature gradient. T L (T R ) denotes the temperature 
maintained at the left (right) lead. Ep — (/xl+//r)/2 de- 
notes the average Fermi energy of the electrodes. Te,r 1 (e) 
and T( tR {e) [T LJ3 = 27r^ k I Ve,/3,k| 2 <5(e - £k)] denote the 
tunneling rates from the QDs to the left and right elec- 
trodes, respectively, e and h denote the electron charge 
and Plank's constant, respectively. For simplicity, these 
tunneling rates are assumed to be energy- and bias- 
independent. Eqs. (1) and (2) have been employed to 
study the thermal properties of single-level QD in the 
Kondo regime. 9 Here, our analysis is devoted to the 
multiple-QD system in the Coulomb blockade regime. 
The expression of the retarded Green function for dot 
£ of a multi-QD system, G\ a (e) can be found in Ref. [7] 
To study the direction-dependent heat current, we 
let T L = T Q + AT/2 and T R = T - AT/2, where 
To = (Tl + T R )/2 is the equilibrium temperature of two 
side electrodes and AT = Tt, — T R is the temperature dif- 
ference. Because the electrochemical potential difference, 
eAV yielded by the thermal gradient could be significant, 
it is important to keep track the shift of the energy level 
of each dot according to eg = Ep + rjgAV/2, where r\i 
is the ratio of the distance between dot I and the mid 
plane of the QD junction to the junction width. Here 
we set i]b = rjc — 0. A functional thermal rectifier re- 
quires a good thermal conductance for AT > 0, but a 
poor thermal conductance for AT < 0. Based on Eqs. 
(1) and (2), the asymmetrical behavior of heat current 
with respect to AT requires not only highly asymmetric 
coupling strengthes between the QDs and the electrodes 
but also strong electron Coulomb interactions between 
dots. To investigate the thermal rectification behavior, 
we have numerically solved Eqs. (1) and (2) for multiple- 
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QD junctions involving two QDs and three QDs for var- 
ious system parameters. We first determine AV by solv- 
ing Eq. (1) with J e — (the open circuit condition) for 
a given AT, T and an initial guess of the average one- 
particle and two-particle occupancy numbers, and eg 
for each QD. Those numbers are then updated accord- 
ing to Eqs. (5) and (6) in Ref. [7] until self-consistency 
is established. For the open circuit, the electrochemical 
potential will be formed due to charge transfer generated 
by the temperature gradient. This electrochemical po- 
tential is known as the Seebeck voltage (Seebeck effect). 
Once AV is solved, we then use Eq. (2) to compute the 
heat current. 

Fig. 2 shows the heat currents, occupation numbers, 
and differential thermal conductance (DTC) for the two- 
QD case, in which the energy levels of dot A and dot B are 
Ea = Ep - AE/5 and E B = E F + a B AE, where a B is 
tuned between and 1. The heat currents are exporessed 
in units of Qo = T 2 /(2h) through out this article. The 
intradot and interdot Coulomb interactions used are Ui = 
30k B T and U ab = T5k B T . The tunneling rates are 
Tar = 0, Y AL = 2I\ and T BR = T BL = T. k B T is 
chosen to be 25r throughout this article. Here, T = 
(r al + Tar)/2 is the average tunneling rate in energy 
units, whose typical values of interest are between 0.1 
and 0.5 meV. The dashed curves are obtained by using 
a simplified expression of Eq. (2) in which we set the 
average two particle occupation in dots A and B to zero 
(resulting from the large intradot Coulomb interactions) 
and taking the limit that T <C k B T so the Lorentzian 
function of resonant channels can be replaced by a delta 
function. We have 

Q/ 1B = n(l-N B )[(l-2N A )(E B -E F )f LR (E B l3) 
+ 2N A (E B + U AB - E F )f LR {E B + U AB )], 

Here Na( B ) is the average occupancy in dot A(B). There- 
fore, it is expected that the curve corresponding to 
E B = E F + 4AE/5 obtained with this delta function 
approximation is in good agreement with the full solu- 
tion, since E B is far away from the Fermi energy level. 
For cases when E B is close to Ep, the approximation is 
not as good, but it still gives qualitatively correct behav- 
ior. Thus, it is convenient to use this simple expression to 
illustrate the thermal rectification behavior. The asym- 
metrical behavior of Na with respect to AT is mainly 
resulted from the condition Tar = and Tal = 2I\ 
The heat current is contributed from the resonant chan- 
nel with e = E B , because the resonant channel with 
e = E B + Uab is too high in energy compared with Ep. 
The sign of Q is determined by fhR{E B )i which indi- 
rectly depends on Coulomb interactions, tunneling rate 
ratio and QD energy levels. The rectification behavior of 
Q is dominated by the factor 1 — 2Na, which explains why 
the energy level of dot-A should be chosen below Ep and 
the presence of interdot Coulomb interactions is crucial. 
The negative sign of Q in the regime of AT < indi- 
cates that the heat current is from the right electrode to 
the left electrode. We define the rectification efficiency as 



t]q = (Q(AT = 30r) - \Q(AT = -30r)|)/Q(AT = 30r). 
We obtain t\q = 0.86 for E B = E F + 2AS/5 and 
0.88 for E B = E F + 4AE/5. Fig. 2(c) shows DTC 
in units of Q k B /T. It is found that the rectification 
behavior is not very sensitive to the variation of E B . 
DTC is roughly linearly proportional to AT in the range 
-20r < k B AT < 20r. In addition, we also find a small 
negative differential thermal conductance (NDTC) for 
E B = Ep + 4A_E/5. Similar behavior was reported in 
the phonon junction system. 10 

Fig. 3 shows the heat current, differential thermal 
conductance and thermal power (S = eAV/k B AT) as 
functions of temperature difference AT for a three-QD 
case for various values of Tar, while keeping T B ( C ).R = 
T B (c),l = T. Here, we adopt r\ A = \Tal - Far|/(2F) 
instead of fixing ija at 0.3 to reflect the correlation of dot 
position with the asymmetric tunneling rates. We as- 
sume that the three QDs are roughly aligned with dot A 
in the middle. The energy levels of dots A, B and C are 
chosen to be Ea = Ep — AE/5, E B — Ep + 2AE/5 
and E c = E F + 3A£/5. U AC = U BA = 15fc B T , 
Ubc = 8k B T 0} Uc — 30k B T , and all other parameters 
are kept the same as in the two-dot case. The thermal 
rectification effect is most pronounced when Tar = 0. 
as seen in Fig. 4(a). (Note that the heat current is 
not very sensitive to U B c)- In this case, we obtain a 
small heat current Q = 0.068Q at AT = -30r, but a 
large heat current Q = 0.33Q at AT = 30r and the 
rectification efficiency t]q is 0.79. However, the heat cur- 
rent for Tar — is small. For Tar = 0.1T, we obtain 
Q = 1.69Q at AT = -30r,Q = 5.69Q at AT = 30r, 
and rjQ = 0.69. We see that the heat current is sup- 
pressed for AT < with decreasing Tar- This im- 
plies that it is important to blockade the heat current 
through dot A to observe the rectification effect. Very 
clear NDTC is observed in Fig. 3(b) for the T A r = 0.1T 
case, while DTC is symmetric with respect to AT for the 
Tar = T A l case. 

From the experimental point of view, it is easier to 
measure the thermal power than the direction-dependent 
heat current. The thermal power as a function of AT is 
shown in Fig. 3(c). All curves except the dash-dotted line 
(which is for the symmetrical tunneling case) show highly 
asymmetrical behavior with respect to AT, yet it is not 
easy at all to judge the efficiency of the rectification effect 
from S for small |AT| (fe B |AT|/r < 10). Thus, it is not 
sufficient to determine whether a single QD can act as an 
efficient thermal rectifier based on results obtained in the 
linear response regime of AT/T < l. 6 According to the 
thermal power values, the electrochemical potential eAV 
can be very large. Consequently, the shift of QD energy 
levels caused by AV is quite important. To illustrate the 
importance of this effect, we plot in Fig. 4 the heat cur- 
rent for various values of Ec for the case with Tar = 0, 
U B c = IO/cbT) and t\a — 0.3. Other parameters are kept 
the same as those for Fig. 3. The solid (dashed) curves 
are obtained by including (excluding) the energy shift 
r/A AV/2. It is seen that the shift of QD energy levels due 
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to AV can lead to significant change in the heat current. 
It is found that NDTC is accompanied with low heat 
current for the case of Eq = Ep + AE/5 [see Fig. 4(b)]. 
Even though the heat currents exhibits rectification ef- 
fect for E c = E F + AE/5 and E c = E F + 3AT/5, the 
thermal powers have very different behaviors. From Figs. 
3(c) and 4(c), we see that the heat current is a highly non- 
linear function of electrochemical potential AV. Conse- 
quently, the rectification effect is not straightforwardly 
related to the thermal power in this system. 

Comparing the heat currents of the three-dot case 
(shown in Figs. 3 and 4) to the two-dot case (shown in 
Fig. 2), we find that the rectification efficiency is about 
the same for both cases, while the magnitude of the heat 
current can be significantly enhanced in the three-dot 
case. For practical applications, we need to estimate the 
magnitude of the heat current density and DTC of the 
IQV junction device in order to see if the effect is signif- 
icant. We envision a thermal rectification device made 
of an array of multiple QDs (e.g. three-QD cells) with a 
2D density N 2 d — 10 n cm~ 2 . For this device, the heat 



current density versus AT is given by Figs. 3 and 4 
with the units Qo replaced by N2dQo, which is approxi- 
mately 9651T/m 2 if we assume T — Q.hmeV . Similarly, 
the units for DTC becomes N 2 dkBQo/T, which is approx- 
imately 34W/°Km 2 . Since the phonon contribution can 
be blocked by the vacuum layer in our design, this de- 
vice should have practical applications near 1A0°K with 
(ksTo « 12.5meV). If we choose a higher tunneling rate 
r > lmeV and Coulomb energy > SOOmeV (possible for 
QDs with diameter less than 1 nm), then it is possible to 
chieve room-temperature operation. 

In summary, we have reported a design of multiplc-QD 
junction which can have significant thermal rectification 
effect. The thermal rectification behavior is sensitive to 
the coupling between the QDs and the electrodes, the 
electron Coulomb interactions and the energy level dif- 
ferences between the dots. 
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Figure Captions 

Fig. 1. Schematic diagram of the isulator/quantum 
dots/vacuum tunnel junction device. 

Fig. 2. (a) Heat current (b) average occupation num- 
ber, and (c) differential thermal conductance as a func- 
tion of AT for various values of Eb for a two-QD junc- 
tion. T A r = 0,T] A = 0.3 and AE = 200r. 

Fig. 3. (a) Heat current, (b) differential thermal con- 
ductance and (c) thermal power as a function of AT for 
various values of Tar for a three-QD junction. 

Fig. 4. (a) Heat current, (b) differential thermal 
conductance and (c) thermal power as functions of AT 
for various values of Ec for a three-QD junction with 
T A r = and r\ A = 0.3. 
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